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摘要 


中 国 十 代 保 留 着 丰富 的 天 文 观 测 记录 ， 其 中 尤为 宝贵 的 是 大 量 实测 的 恒星 观测 数据 。 然 而 ， 确 定 这 些 观测 数据 的 
年 代 是 一 个 相当 复杂 的 问题 。 本 文 以 宋 元 时 代 的 两 组 观测 数据 为 基础 ， 运 用 广义 霍 夫 变换 的 技术 手段 对 数据 进行 
计算 和 分 析 。 我 们 成 功 地 获得 了 对 恒星 观测 年 代 的 较为 准确 的 估算 结果 。 这 一 方法 为 未 来 分 析 更 多 古代 天 文 星 表 


数据 提供 了 新 的 途径 和 视角 。 
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1. 简介 


中 国 古 代 有 丰富 的 天 文 观测 记录 ， 尤 其 珍 贯 的 是 有 
不 少 实测 的 恒星 观测 记录 。 近 代 以 来 ， 不 断 有 学 者 对 这 
些 数据 进行 分 析 ， 以 确定 这 些 数据 的 观测 年 代 等 。 传 统 
的 方法 有 岁差 改正 法 册 1， 利用 岁差 的 数据 对 这 些 观测 数 
据 进 行 逼近 分 析 。 还 有 傅立叶 分 析 方 法 趾 ， 利 用 恒星 数 
据 对 赤 纬 误差 进行 时 代 的 负 算 。 本 文采 用 了 广义 霍 夫 变 
换 对 古代 恒星 观测 数据 进行 了 处 理 ， 并 获得 了 较 好 的 年 
代 佑 计 。 

之 前 学 者 们 确定 年 代 的 方法 需要 有 较为 准确 的 准 
确 的 古今 恒星 对 应 数据 和 恒星 的 位 置 坐标 ， 本 文 将 利用 
广义 霍 夫 变换 的 原理 来 实现 一 个 新 的 算法 来 提升 计算 
的 精度 ， 广 义 霍 夫 变换 使 用 的 的 是 概率 计算 ， 不 需要 完 
备 准确 的 数据 ， 仅 需要 大 部 分 数据 相对 正确 即 可 ， 在 计 
算 中 也 仅 需要 恒星 的 去 极度 数据 ， 而 其 他 方法 不 仅 需 要 
去 极度 数据 也 需要 和 人 入 度数 据 。 

本 文 将 探讨 使 用 广义 霍 夫 变换 的 技术 方法 来 确定 
中 国 古代 恒星 的 观测 年 代 。 


2. 技术 方法 


霍 夫 变换 (Hough Transform ) 是 一 种 特征 提取 技 
术 ， 被 广泛 应 用 于 图 像 分 析 、 计 算 机 视觉 领域 。 最 初 是 
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由 霍 夫 《Hough) 于 1962 年 提出 的 。 它 的 原理 基于 以 下 
思想 : 对 于 图 像 中 的 每 个 点 ， 其 在 参数 空间 中 可 能 是 某 
种 形状 ， 但 当 所 有 的 点 都 描述 在 参数 空间 时 ， 通 过 类 似 
于 投票 的 机 制 就 可 以 求 出 所 需 的 参数 。 也 就 是 说 ， 霍 夫 
变换 是 一 个 概率 算法 ， 通 过 参数 空间 的 分 布 概率 计算 来 
确定 参数 。 广 义 霍 夫 变换 是 传统 霍 夫 变换 的 一 个 扩展 。 

对 于 本 文 的 问题 求解 来 涪 ， 恒 星 观测 时 期 的 北极 点 
的 位 置 就 是 所 求 的 参数 ， 位 置 坐标 数据 组 成 的 就 是 参数 


空间 。 


2.1. 算法 原理 


图 1: 恒星 坐标 示意 图 


中 国 古 代 的 天 球 坐 标 系 统 是 世界 上 最 早 的 亦 道 坐 
标 系 统 ， 天 文学 家 在 进行 恒星 观测 时 ， 和 现代 的 天 文 观 
测 类 似 ， 首 先 需要 定好 北极 。 与 现代 恒星 的 位 置 数据 类 
似 ， 中 国 古 代 恒 星 观测 数据 也 由 两 部 分 组 成 : 入 箱 度 和 
去 极度 。 

在 传世 文献 中 对 恒星 观测 数据 的 标注 格式 一 般 为 
“入 某 宿 若干 度 若干 分 ”和 “去 极 若干 度 若干 分 "。 以 
图 1 为 例 ，X 点 是 某 答 距 星 的 位 置 ，P 是 北 天 极 ， 那 么 
PX 标注 着 这 一 答 的 基准 经 线 ， 对 于 该 答 的 一 颗 恒 星 9 
来 说 ， 经 度 差 9 就 是 人 宿 度 ， 北 天 极 到 S 点 的 角度 w 
就 是 去 极度 ， 其 中 去 极度 与 现代 赤道 坐标 邓 的 杰 纬 关系 
为 : 


6=90—w 


作为 对 比 ， 现 代 天 文 恒 星星 表 中 包含 的 天 体位 置信 
上 息 ， 一 般 使 用 赤道 坐标 系 ， 定 义 为 赤 经 和 杰 纬 : 


s = (a,,6,) 


去 极度 数据 w 是 恒星 与 北 天 极点 的 角 距 离 。 恒 星 
观测 时 北 天 极点 的 位 置 设 为 P = (a,6)。 那 么 对 任 一 颗 
恒星 与 并 的 当时 的 北极 点 来 说 ， 向 量 s 和 p 的 点 积 为 : 


s.:p=|sllplcosw 四 
一 Cos 

在 需要 求解 的 北极 点 (a, 6) 参数 空间 中 ， 每 组 数据 
对 应 的 北极 点 ， 在 参数 空间 可 以 表示 为 一 个 圆 ， 如 图 2， 
这 里 选取 了 9 个 北极 点 附近 的 较 亮 恒星 观测 数据 作为 
演示 ， 可 以 清晰 的 看 到 参数 空间 的 圆 在 一 个 点 附近 进行 
有 汇聚 ， 只 有 一 颗 星 稍微 远离 了 汇聚 点 。 当 有 更 多 的 数 
据 时 ， 可 以 根据 概率 获得 最 佳 的 数据 汇聚 点 ， 也 就 是 北 
极点 的 候选 位 置 。 

那么 就 可 以 通过 广义 霍 夫 变换 的 计算 原理 ， 对 所 有 
恒星 观测 数据 的 处 理 计算 可 以 最 终 获 得 可 能 的 北极 点 
坐标 , 然后 通过 北极 点 的 位 置 以 及 岁差 数据 (如 图 3) 就 
可 以 确定 恒星 观测 的 确切 年 代 。 

由 于 地 球 的 岁差 ， 地 球 的 自转 轴 有 个 大 约 26000 年 
的 转动 周期 。 公 元 前 500 年 至 2100 年 的 北极 点 如 图 3 所 
示 。 本 文 计算 北极 点 历史 的 位 置 使 用 了 Jan Vondr&k 中 
的 岁差 模型 算法 。 


图 2: 霍 夫 变换 参数 空间 


Precession (-500 to 2100) 


22h 85° 23h Oh 1lh 2h3h4h5h 6h A 85° 8h 


图 3: 公元 前 500 年 至 2100 年 北极 点 的 位 置 


2.2. 算法 步 又 


在 这 里 假定 一 个 数据 集 的 观测 数据 是 在 一 个 较 
短 的 时 期 进行 的 ， 这 种 情况 下 ， 可 以 设 定 北极 点 的 
移动 可 以 忽略 。 这 组 观测 数据 在 J2000 历 元 下 的 的 
坐标 为 (qa,,6)(i = 1,2,...,m)， 转 换 为 第 卡尔 坐标 为 
(Xj, yi,2Z0)(i = 1,2,...,m)。 它 们 的 去 极度 数据 为 wj(i = 
1,2,...,m)o 

假设 这 一 时 期 的 北极 点 位 置 为 (w,56)， 转 换 为 笛 卡 
尔 华 标 为 p(x, y, z)。 

任 取 两 颗 恒 星 A 和 B 的 数据 ,空间 位 置 (xj 2) 
和 sj(xj,yj,2zj)， 基 中 i 关 j。 那 么 有 


sj*P=cosw®; (2) 
Ip|=1 


在 笛 卡 尔 坐 标 系 下 : 


XX + yy+ Zz = COS Qi 
XjX+yy+2ZiZ= cos; (3) 
x +y+2z =1 
求解 上 述 方 程 组 ， 如 图 4， 所 示 ， 这 个 方程 组 的 解 
有 三 种 情况 : 唯一 和 解 、 双 解 和 无 解 。 那 么 通过 和 迭代 ， 可 
以 计算 出 任 两 个 恒星 之 间 的 交点 解 。 
由 于 观测 数据 有 一 定 的 误差 ， 在 求解 方程 组 的 过 程 
中 ， 可 以 进行 蒙特 卡 洛 随机 抽样 模拟 ， 可 以 有 效 地 增加 
结果 的 精确 度 。 


图 4: 方程 组 求解 


2.3. 原始 数据 预 处 理 


在 数据 处 理 分 析 前 ， 还 需要 对 去 极度 数据 进行 换 
算 。 观 测 数据 中 的 度数 是 中 国 传统 的 古 度 系统 ， 即 一 个 
周 天 365.25 度 ， 因 此 在 做 数据 处 理 时 还 需要 将 数据 转 
换 为 周 天 360 度 的 角度 系统 。 


360 


1 
十 度 365.25 


degree 


3. 参考 数据 


精确 的 结果 来 源 于 高 精度 的 星 表 ， 在 数据 处 理 过 程 
中 ， 我 们 用 用 了 最 新 的 天 体 测 量 框架 和 最 新 的 恒星 观 
测 数据 作为 数据 计算 的 基础 。 本 文采 用 的 计算 库 为 IAU 
的 基础 天 文学 标准 库 (SOFADP' 9 )。 

现代 恒星 星 表 数据 采用 了 亮 星 数 据 比 较 精 确 的 
1997 年 发 布 的 伊 巴 谷 星 表 []， 然 后 为 参考 了 2007 年 
新 处 理发 布 的 伊 巴 谷 星 表 凶 以 及 2012 年 以 依 巴 谷 星 表 
为 基础 的 扩展 伊 巴 谷 星 表 XHIPM。 后 者 除了 新 处 理 之 
外 ， 还 加 入 了 最 新 的 恒星 自行 、 视 差 和 视 向 速度 等 数 
据 ， 通 过 合并 这 几 个 数据 源 ， 我 们 构建 了 一 个 服务 于 中 
国 古 代 星 表 计 算 的 基准 星 表 : 中 国 古 代 天 文 基础 参考 星 
表 00， 星 表 的 数据 历 元 采用 J2000.0。 基 于 这 些 数据 和 
岁差 计算 框架 ， 可 以 相对 精确 地 计算 不 同 历 史 时 期 恒星 
的 位 置 数据 。 

算法 步骤 是 初步 估算 一 个 古代 星 表 观测 的 年 代 后 ， 
再 根据 自行 等 数据 再 迭代 计算 星 表 中 恒星 在 观测 年 代 
的 位 置 ， 以 最 终 确认 数据 的 观测 年 代 。 


4. 数据 处 理 和 分 析 


本 文 以 宋代 和 元 代 两 份 星 表 为 例 对 上 述 算法 进行 
测试 和 分 析 。 


4.1. 宋代 星 表 

北宁 年间， 社会 经 济 和 科学 文化 高 度 发 达 ， 在 天 文 
历法 方面 ， 多 次 制作 核 正 天 文 仪器 ， 并 对 星象 进行 了 大 
量 的 观测 ， 其 成 果 和 记录 在 历代 都 是 罕见 的 ， 是 中 国 古 
代 科 技 的 一 个 高 峰 时 期 。 在 现存 传世 文献 中 ，“ 灵 吝 秘 
苑 >、<《 文 献 通 考 >、《“《 管 宕 辑 要 >、《 陪 象 通 览 》 和 《天 元 
历 理 > 等 都 保存 有 大 量 的 恒星 观测 记录 。 潘 得 等 U1 站 
鲁 对 中 国 国 家 图 书馆 藏 明 抄 本 “《“ 灵 人 训 秘 苑 > 以 及 其 他 的 
传世 文献 进行 整理 分 析 ， 编 辑 出 了 包含 360 颗 星 的 恒星 


星 表 ， 并 定名 为 《 宋 皇 佑 星 表 >。 本 节 以 该 星 表 为 基础 ， 
以 广义 霍 夫 变换 的 算法 对 其 再 进行 处 理 ， 必 重新 定 其 年 
代 。 

计算 结果 如 表 1 所 示 ， 其 中 极 益 是 候选 北极 点 与 最 
近年 份 北极 点 的 为 角度 差 。 在 图 5 中 ， 可 以 看 到 选取 不 
同 的 误差 范围 ， 可 以 得 到 不 同 但 近似 的 候选 结果 ， 其 中 ， 
o = 0.25 和 o = 0.3 上 古 度 的 情况 下 ， 候 选 点 收敛 集中 的 
程度 较 好 ， 图 6 是 o = 0.25 情况 下 的 图 形 。 因 此 可 以 确 
定 的 年 代 大 约 为 公元 1011 至 1016 年 前 后 。 
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图 5: 计算 结果 的 散 点 图 


公元 1011 年 为 大 中 冬 符 四 年 ， 据 《灯会 要 > 中 记 
载 :“ 真 宗 大 中 祥 符 三 年 间 二 月 四 上 日 ， 司 天 监 言 : ' 冬 官 
正 韩 显 符 造 铜 浑 仪 成 。 诏 移入 龙 图 阁 ， 令 显 符 选 学 生 
中 可 教 者 传授 其 业 ”。 同 书 还 记载 有 :“' 祥 符 三 年 七 月 五 
日 ， 司 天 监 冬 官 正 韩 显 符 上 外 官 星 位 去 斗 极 度数 。 棠 文 
目 宿 曜 度 分 域名 录 一 卷 2?2。 这 一 信息 同时 被 摘录 入 “ 玉 
海 > 0 “ 衬 符 星 度 ”条 目下 。 

另 据 史料 记载 ， 韩 显 符 还 参加 了 至 道 元 年 〈 公 元 
995 年 )“ 至 道 司 天 宫 铜 浑 仪 ” 的 制作 ， 他 在 制作 浑 仪 的 
同时 ， 还 留 下 了 详细 的 制作 方法 ，< 宋 史 . 天 文 志 > [1 
记载 有 :“ 铜 候 仪 ， 司 天 冬 官 正 韩 显 符 所 造 ， 其 要 本 淳 风 
及 僧 一 行 之 遗 法 。 显 符 自 著 经 十 卷 上 之 书 府 。 铜 仪 之 制 
有 九 ”。 另 韩 显 符 于 大 中 冬 符 六 年 去 世 ， 共 学 生 宋 代 天 
文学 家 杨 惟 德 代 留 泰 名 编 扎 “ 景 佑 喜 象 新 书 >， 在 编写 


此 书 的 过 程 中 ， 杨 惟 德 等 人 也 全 
测 。 


进行 过 一 次 周 天 恒 星 观 


o=0.250 Year:1011 Diff:0.04° 


图 6: o = 0.25 (〈 古 度 ) 极点 位 置 


至 此 大 致 可 以 推测 公元 1010 年 ， 冬 官 正 韩 显 符 监 
造 了 冬 符 龙 图 图 铜 浑 天 仪 ， 甘 守 分 其 教授 观测 ， 并 留 下 
了 观测 记录 一 涂 。 基 学生 杨 惟 德 继续 韩 显 符 的 事业 ， 也 
组 织 过 恒星 观测 。 传 世 文 献 中 记录 的 数据 有 可 能 是 此 次 
观测 的 记录 ， 并 在 随后 的 几 十 年 间 的 若干 次 观测 中 被 部 
分 的 承袭 了 下 来 。 


4.2. 元 代 星 表 

中 国 国 家 图 书馆 藏 的 一 部 明 抄本 “天 文 莫 钞 >》09， 
共 十 一 种 六 册 ， 其 中 第 三 册 为 《三 垣 列 舍 人 往 去 极 集 > 
一 餐 ， 记 录 了 约 740 颗 恒 星 的 详细 观测 记录 。 上 世纪 八 
十 年 代 被 整理 发 现 后 ， 多 名 学 者 对 此 进行 了 研究 。 陈 
应 上] 在 1986 年 对 该 数据 进行 了 分 析 ， 结 诊 是 该 星 表 观 
测 于 元 代 公 元 1280 年 前 后 ， 那 个 时 代 恰 好 是 郭守敬 建 
造 浑 仪 的 时 代 ， 因 此 ， 他 认为 这 是 郭守敬 进行 恒星 观测 
的 结果 。 潘 才 口 也 全 研究 分 析 过 该 数据 ， 结 诊 与 陈 文 类 
似 ， 御 定年 代为 公元 1280 年 左右 。 

上 世纪 九 十 年 代 孙 小 沾 启 也 对 该 星 表 进行 了 严密 
的 数据 分 析 后 ， 测 定 其 实际 观测 年 代为 明 初 洪武 年 间 的 
公元 1380 年 ， 并 发 表 了 其 中 678 颗 恒 星 的 证 认 结 果 。 

本 文 以 孙 小 淳 文中 发 表 的 669 颗 认 证 恒星 数据 〈 去 
掉 了 部 分 不 完整 数据 ) 为 基础 进行 计算 ， 去 极度 的 误差 


表 1: 宋代 星 表 处 理 结果 


数据 源 核定 值 灵 宇 秘 苑 一 。 灵 秘 苑 二 文献 通 考 象 林 管 宪 辑 要 天 元 历 理 
星 数 352 342 347 263 297 340 300 
去 极度 误差 ”年代 ” 极 差 年 代 极 差 年 代 极 差 年 代 极 差 年 代 极 差 年 代 极 差 年 代 极 差 
( 古 度 ) (deg) (deg) (deg) (deg) (deg) (deg) (deg) 
0.10 995 0.03 989 0.11 987 0.07 1003 0.07 998 0.07 983 0.11 1002 0.08 
0.20 1005 0.04 996 0.01 996 0.06 1015 0.06 1006 0.05 997 0.15 1012 0.07 
0.25 1011 0.04 1002 0.00 1002 0.06 1022 0.05 1012 0.04 1009 0.16 1019 0.06 
0.30 1016 0.05 1007 0.01 1008 0.07 1028 0.04 1018 0.03 1019 0.17 1026 0.05 
0.35 1019 C0.05 1010 0.02 1012 0.07 1033 0.04 1025 0.03 1026 0.18 1033 0.04 
0.40 1021 0.05 1013 0.03 1016 0.07 1034 0.03 1028 0.03 1029 0.17 1036 0.04 
0.45 1022 0.06 1013 0.03 1019 0.08 1033 0.02 1032 0.03 1030 0.16 1036 0.03 
0.50 1022 0.06 1014 0.04 1020 0.08 1035 0.02 1031 0.03 1030 0.16 1036 0.03 
范围 以 高 斯 分 布 进行 蒙特 卡 洛 抽样 ， 计 算 的 结果 如 表 2。 
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图 7: 计算 结果 的 散 点 图 图 8: o = 0.1( 古 度 ) 极点 位 置 
通过 表 2 和 图 8、 图 7， 我 们 可 以 看 到 ， 在 o = 0.1 上古 
度 的 情况 下 ， 数 据点 收敛 的 比较 好 。 可 以 确定 的 年 代 大 
约 为 公元 1355 年 前 后 。 表 2: 《天 文 汇 抄 》 星 表 
公元 1355 年 是 元 朝 的 至 正 十 五 年 ， 这 已 经 是 元 朝 去 极度 误差 十 度 ) 北极 年 代 极 基 (deg) 
， 、 、 0.050 1353 0.20 
的 末期 了 。 从 数据 的 记录 来 看 ， 这 是 一 份 及 其 难得 的 数 ee i 
据 记录 ， 数 据 以 星 图 形式 记录 ， 在 恒星 侧 边 记 录 基 位置 0.125 1357 0.20 
y y > 、 a » 0.200 1361 0.20 
数据 。 另 外 其 数据 的 分 度 为 百分制 ， 这 在 之 前 也 是 没有 
0.250 1363 0.20 
0.333 1366 0.20 


过 的 。 可 以 推测 ， 这 应 该 是 元 代 末 期 的 司 天 监 官员 使 用 
郭守敬 所 造 浑 仪 进行 的 一 次 细致 的 恒星 观测 ， 并 且 以 图 
的 形式 的 保存 下 来 ， 而 且 在 明代 有 一 次 传 抄 。 


5. 小 结 和 讨论 


通过 上 壕 两 组 数据 的 计算 和 分 析 ， 我 们 可 以 初步 对 
个 广义 霍 夫 变换 以 确定 年 代 的 算法 进行 评估 ， 该 算法 可 
以 有 效 的 估算 中 国 古 代 恒 星 的 观测 年 代 ， 获 得 了 上 比 之 前 
研究 更 好 的 精确 度 。 和 其 他 确定 年 代 的 方法 相 比 ， 使 用 
广义 霍 夫 变 换 的 优点 是 可 以 不 用 考虑 赤 经 数据 ， 也 不 用 
特别 考 虚 错 误 的 数据 ， 在 大 规模 的 误差 条 样 后 ， 可 以 有 
效 消除 错误 数据 的 影响 。 同 时 在 计算 过 程 中 ， 也 充分 考 
虑 了 恒星 的 自行 等 相关 参数 ， 进 一 步 提 高 了 年 代 确定 的 
精度 。 

在 本 文 的 计算 过 程 中 ， 灯 代 观 测 数 据 需 要 在 c = 
0.25 的 情况 下 ， 数 据 聚 集 度 较 好 ， 而 到 了 元 代 ，c = 0.1 
的 情况 下 ， 数 据 已 经 聚集 度 就 超过 了 宋代 。 这 显示 出 从 
宋代 到 元 代 ， 随 着 仪器 的 进步 ， 数 据 精 度 也 在 明显 的 提 
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本 文 使 用 的 方法 在 宋 元 星 表 上 取得 了 较 好 的 计算 
结果 ， 这 也 会 为 下 一 步 分 析 更 多 的 古代 天 文 星 表 数据 提 
高 了 一 个 新 的 方法 和 思路 。 
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